Crossovers between superconducting symmetry classes 
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We study the average density of states in a small metallic grain coupled to two superconductors 
with the phase difference n, in a magnetic field. The spectrum of the low-energy excitations in the 
grain is described by the random matrix theory whose symmetry depends on the magnetic field 
strength and coupling to the superconductors. In the limiting cases, a pure superconducting sym- 
metry class is realized. For intermediate magnetic fields or couplings to the superconductors, the 
system experiences a crossover between different symmetry classes. With the help of the supersym- 
metric er-model we derive the exact expressions for the average density of states in the crossovers 
between the symmetry classes A-C and CI-C. 

PACS numbers: 74.45.+C, 73.22.Dj, 73.22.Gk 



Introduction. Energy levels in small metallic parti- 
cles with chaotic electron dynamics are random numbers. 
It is generally accepted that their spectral statistics in 
the ergodic regime is described by the random matrix 
theory (RMT) [lj. For disordered grains, this had been 
proved by Efetov Q with the help of the supersymmetry 
technique while for quantum billiards in the absence 
of disorder this statement is usually referred to as the 
Bohigas conjecture 0. 

In the RMT, a system is characterized solely by its 
symmetry. In the application to condensed matter, 
the standard three Wigner-Dyson ensembles (orthogo- 
nal, unitary and symplectic) |5( describe level statistics 
in small metallic grains in the presence or absence of the 
time-reversal and spin-rotation symmetries . 

Recently, the Wigner-Dyson classification had been ex- 
tended to superconducting and chiral [3j symmetry 
classes, which arise when the Hamiltonian possesses an 
additional symmetry with respect to changing the sign of 
the energy (counted from the Fermi energy). With the 
appearance of a selected energy point, in the supercon- 
ducting/chiral classes even the average density of states 
(DOS), (p(E)), becomes a nontrivial function of the en- 
ergy. This should be contrasted to the standard Wigner- 
Dyson ensembles where (p{E)) = S^ 1 = const and the 
first nontrivial quantity is the pair correlation function 
Ra(u) = S 2 (p(E + u)p{E)) - 1. 

The symmetry classes (three Wigner-Dyson, four su- 
perconducting and three chiral) correspond to the limits 
when various symmetries are either present or completely 
broken. In the intermediate cases, the system experiences 
a crossover between different symmetry classes. The pair 
correlation function in the crossover between the orthog- 
onal and unitary classes was obtained in Refs. [H, 
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FIG. 1: A normal-metal dot coupled to two superconducting 
terminals with the phase difference ir, in a magnetic field. NS 
interfaces are characterized by the set of transparencies {Ti}. 



where x = ttuj / S and a is the symmetry-breaking param- 
eter. Equation (JTJ interpolates between the orthogonal 
(a = 0) and unitary (a — oo) results. The pair corre- 
lator in known also in the symplectic-unitary crossover 
fiol [ll| , its form being similar to Eq. (JTJ . 

The purpose of this Letter is to theoretically study 
crossovers between superconducting classes. 

We will calculate the average density of states in a 
small diffusive metallic grain coupled to two supercon- 
ducting terminals through tunnel barriers, see Fig. [T] 
The terminals have the phase difference 7t ensuring the 
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FIG. 2: Crossovers between spin-symmetric symmetry classes 
driven by the magnetic field (H) and coupling to supercon- 
ductors (A). The dimensions of the FF and BB sectors of the 
supersymmetric <r-model for the average density of states are 
shown by (nF,ns). 
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absence of the minigap in the excitation spectrum [12j . 
A magnetic field H is applied to the system. The spin- 
rotation symmetry is assumed to be intact. We will 
be interested in the ergodic regime, E <C £rh> where 
-Erh = D/L 2 is the Thouless energy, D is the diffusion 
constant, and L is the grain size. 

Under these conditions, the excitation spectrum in 
the grain can be described in terms of the RMT in 
the crossover region between the four symmetry classes 
shown in Fig. [2j 

Mapping to Efetov's cr-model. First attempts of 
field-theoretical description of hybrid NS systems 13l-ll7 
inspired by the identification of superconducting sym- 
metry classes @ have used the Bogolyubov-de Gennes 
(BdG) Hamiltonian as the starting point, 



TABLE I: Basic matrices in the two versions of the cr-model. 
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where H is the single-particle Hamiltonian, and A(r) is 
the pairing field. The average quasiparticle DOS, 



(p(£,r)} = -lm<trel(r))/7r, 



(3) 



is expressed in terms of the retarded Green function of 
the BdG Hamiltonian, 



g% = (E-n BdG + zoy 



(4) 



which is then represented as a functional integral over an 
8x8 supermatrix field Q acting in the direct product 
FB <g> N ® PH of the Fermi-Bose (FB), Nambu (N) and 
Particle-Hole (PH) spaces (spin-symmetric case is con- 
sidered) . 

In hybrid NS systems, Andreev reflection off the order 
parameter field A couples the states with opposite en- 
ergies, E and —E. So, the Nambu-Gor'kov Green func- 
tion Ge essentially involves a pair of the retarded and 
advanced normal-metal Green functions, and G^ E . 
In the absence of the superconducting pairing field, A(r), 
correlations between the latter are conveniently described 
by Efetov's supersymmetric cr-model 0, Q of the orthog- 
onal symmetry class, with an 8 x 8 superfield Q acting 
in the direct product FB (g) RA <g> TR of the Fermi-Bose 
(FB), Retarded- Advanced (RA) and Time-Reversal (TR) 
spaces (again we assume no spin interactions). 

Thus, in studying the proximity effect in the normal 
part of a hybrid system, it is tempting to reformulate 
the field theory of Refs. [IMll 

in the language of Efe- 
tov's supersymmetric cr-model. Provided that the inverse 
proximity effect in the superconducting regions can be 
neglected (rigid boundary conditions), Andreev scatter- 
ing of normal electrons off the superconducting terminal 
will be viewed as an effective boundary condition at the 
NS interface mixing the R and A components of the field 
Q. Such a description is close in spirit to the scattering 
approach [18j . 
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The average local DOS is given by the functional inte- 
gral over the normal- metal region 171 ]: 



{p{E,r)) = ^Re^ sti(k^Q)e~ SDlQ] ~ SrlQ] DQ(r), (5) 

where Sd[Q] is the bulk action: 

s D = y J dr str [ D ( V< 3 + * eA [ r 3,Q]) 2 + 4z£;AQ], 

(6) 

and the action Sr[Q] = Sri[Q] + SValQ] describes NS 
interfaces (il. [l9|: 

S r« = -^E^Ml + e-^Q^Q^]- (7) 



Here v is the DOS per one spin projection at the Fermi 
level, A is the vector potential, Q( a > labels the Q field 
at the boundary with the a-th superconductor, and the 
NS interface is specified by transmission coefficients Tj = 
l/cosh 2 /^, with i labelling open channels. The field Q 
satisfies Q 2 = 1 and is subject to an additional symmetry 
constraint 



Q = CQ 1 C 



(8) 



In the NS cr-model for (£/f ) [ij}, the matrices A, r 3 
and C are given by the first column of Table HI An exact 
mapping to Efetov's cr-model is realized by the similarity 
transformation Q i— > VQV^ 1 with the matrix 
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(9) 



where the inner (outer) grading corresponds to the PH 
(N) space, and k — diag(l, — 1)fb (we follow notations of 
Ref. |3|). Conjugation by V simultaneously transforms 
the matrices A, T3 and C from NS representation to Efe- 
tov's representation given by the last column of Table |U 
This provides an exact mapping between the NS cr-model 
for Green function of the BdG Hamiltonian, (Ge), to 
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the standard Efetov's orthogonal cr-model for the prod- 
uct (G^G^ E ). On such a mapping both the structure 
of the manifold and the cr-model action get reproduced. 
We emphasize that this mapping takes place only in the 
normal part of a hybrid NS system, where the pairing 
amplitude A = 0. 

To complete the formulation of the model we have to 
specify the Q matrix in the bulk of a superconductor, 
Qs- It has a familiar form parameterized with the help 
of the spectral angle 6s = arctan(iA/£') as 

Qs = A cos 9 s + t sin 6 S . (10) 

The most nontrivial ingredient of the mapping from the 
SN cr-model to Efetov's cr-model is the form of the matrix 
S. In the initial NS representation [17| it is just the Pauli 
matrix in the Nambu space: cr^. Conjugating by V we 
get it in Efetov's representation: 

£ =(A o) • Hit) ' (n) 

\ / TR V / RA 

We see that superconducting boundary conditions violate 
super symmetry: The matrix £ acts as af A aJ R in the FF 
block and as cripV™ in the BB block. This is the reason 
why a nontrivial DOS can be obtained by integration ([5]) 
over the standard orthogonal cr-model manifold. 

Zero-dimensional limit. In the ergodic regime, E <C 
i?Th, the functional integral (|SJ) is dominated by the zero 
mode, Q(r) = const. We will be interested in the average 
global DOS normalized by the inverse mean quasiparticle 
level spacing 8 = (2vV)~ 1 : 

( e (E))=sJ{p(E,r))dr. (12) 

This quantity can be written as an integral over a single 
8x8 supcrmatrix Q: 

(g(E)) = 1 Re J str (k AQ) e~ s [Q] DQ, (13) 

with the action consisting of three terms: 

S[Q] = — str AQ - - str(r 3 Q) 2 + ^ str(SQ) 2 . (14) 
4 4 8 

Here x = nE/8, and the symmetry breaking parameters 
a and 7 are given by 

a = iivDe 2 J A 2 c?r, 7 = (15) 

One can estimate a ~ (<f> / '4>o) 2 1 9 , where c6 is the flux 
through the grain, 0o is the flux quantum, and g ~ 
Etyi/S ^ 1 is the dimensionless grain conductance. The 
last term in the action (fM)) is written in the tunneling 
limit, Tj <C 1, and the parameter 7 is expressed through 
the dimensionless (in units of e 2 /h) Andreev conductance 



of the grain [20] , G a = 2 J2i (the factor 2 accounts for 
two NS interfaces assumed to be identical, see Fig. [TJ. A 
strong magnetic field randomizes electron phase and the 
crossover from the unitary class to class C can be ob- 
tained with just one superconducting terminal attached 
to the grain (in that case Ga = The last term 

in Eq. (fT4|) is written in the subgap limit, E <C A, when 

Equations (fT3")) and describe the average DOS in 
the two-parametric crossover between the four symmetry 
classes shown in Fig. [3] Instead of studying the compli- 
cated general behavior, we will restrict ourselves to two 
one-parametric crossovers: class A - class C (a — 00, 7 
arbitrary) and class CI - class C (7 = 00, a arbitrary). 

Class A — class C crossover. In sufficiently strong 
magnetic fields (a 3> 1), cooperon degrees of freedom 
get frozen and the Q matrix becomes diagonal in the TR 
space: 

Q =( Q 0*MOV) TR ' ™ 

where the 4x4 matrix Q u G FB ® RA spans the mani- 
fold of the unitary Efetov's cr-model. Therefore one can 
simply take the well-known Efetov's parametrization of 
this manifold in which the FF and BB sectors are 
parametrized by a radial variable (Af,b) and an angular 
variable (c/?f,b), with —1 < Ap < 1 and Ab > 1. It can 
be easily seen that the term (7/8) str(SQ) 2 in the action 
(fl4]) does not depend on Af. Thus coupling to a super- 
conductor suppresses only the BB degrees of freedom, 
shrinking it to a point at 7 — > 00 (see Fig. [2|) . Calcula- 
tion of the integral (fl~3|) is straightforward leading to the 
exact expression for the class A - class C crossover: 

/>oo 

= 1 - 27— / dXX cos Ax e-^ 2 - 1 ). (17) 




FIG. 3: The average DOS, (g(E)), in the class A (unitary) - 
class C crossover. The curves correspond to different values of 
the symmetry-breaking parameter: 7 = (unperturbed DOS, 
(g{E)} = 1), 0.01, 0.1, 1, and 00 (class C). 
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The function (g(E)) is plotted in Fig.[3]for several values 
of the symmetry breaking parameter 7. 

In the limit 7 — > 00, one recovers the C-class result: 



(qc(E)) = 1 - 



sin 2x 
2x 



(18) 



In the limit of weak coupling to a superconductor, 7 <C 
1, the uniform metallic average DOS is perturbed in a 
small vicinity of the Fermi energy, at E < <5^/7, but this 
perturbation is strong, completely suppressing the DOS 
at E = 0: 



(9(E)) « / 



2^7 



/(*) - 2ze- 



(19) 



Vanishing of (g(0)) can be explained by repulsion of en- 
ergy levels E and —E which becomes effective at very 
small energies E < 8 yFy. (Formally, noncommutativity 
of the limits E — > and 7 — > is a consequence of the 
noncompactness of the BB sector of the theory.) 

Level statistics in classes A and C is known to be de- 
scribed by free fermions: The joint probability density 
of energy levels can be interpreted as the square of the 
ground-state wave function for a system of nonintcract- 
ing one-dimensional (ID) fermions. Class A (unitary) 
corresponds to free fermions on a line l], while class C 
corresponds to free fermions with the Dirichlet bound- 
ary condition at the origin Q. Our result (TTTl) indicates 
that in the crossover A-C energy levels cannot be consid- 
ered as noninteracting fermions in an appropriate single- 
particle potential. 

Class CI — class C crossover. Now we turn to 
the case of strong coupling to a superconductor, 7 ^> 
1. The corresponding term in the action (fl4| enforces 
str(EQ) 2 = 0. The latter condition together with the 
relation {A, £} = leads to the linear constraint 



{Q,£} = 0. 



(20) 



The next step in constructing the parametrization suit- 
able for calculation in the crossover region is to study the 
commuting generators W of the Cl-class manifold. Writ- 
ing Q = A(l + W + ... ) with {W, A} = 0, and solving 
the linearized constraints (jSJ and (|20[) , we find four com- 
muting generators, three residing in the FF sector, and 
one in the BB sector (see Fig. [2]): 



(w FF 



V 



' f com 1 ' r com' 



(21) 



where 



W FF = - 



( z -c\ 
c z* 
z* c 

V-c z J 



(22) 



W B B = ~ 



/ooo l\ 

10 
10 
\i 0/ 



(23) 



RA 



with complex z and real c and q. 

As the a-dependent term in the action (|14[) contains 
the symmetry-breaking matrix T3, it is convenient to split 
W com into a part Wcom commuting with T3 (z and z* 
modes) , and a part Wcom anticommuting with T3 (c and 
q modes). At finite a, the latter modes acquire a mass 
proportional to a. They will completely freeze out in the 
C-class limit (a — > 00), where only the modes z = 9e lip 
and z* = 6e~ lv will be unaffected, generating the sphere 
S 2 in the FF sector. 

In constructing the global parametrization of the CI- 
class manifold we will follow an approach of Ref. [9J in 
order to maximally simplify the symmetry breaking term 
S a [Q] — — (a/4) str(r 3 <5) 2 . We start with parametrizing 
the commuting content of Q as Q CO m = Uq 1 U~ 1 AU rn Uc, 
where 



U c = exp W«, U m = exp W ( 



(-) 



(24) 



With such a choice, S a [Qcom] will explicitly depend only 
on the massive modes c and q. 

Now we turn to the Grassmann content of the 
parametrization. We search for Grassmann generators 
W which obey [T3, W] = (to simplify the term S a ) and 
[A, W] = (to simplify the term strAQ). Employing 
also linearized Eqs. (JSJ) and (|2"0"|) . we find that these an- 
ticommuting generator parametrized by two Grassmann 
numbers: 




FB 




(25) 



(26) 



FB 



The desired parametrization of the Q manifold in the 
spirit of Ref. [9J has the form 



Q 



U^Uj-U^U^KUJUpUcUt, 



(27) 



where the matrices U m and Uc are defined in Eq. 
and 



U t = exp W G + J K, p], U,= exp W£> [fi, V ] . (28) 

After some algebra we obtain the Berezinian of the 
parametrization (I27|) : 



(+)[ 



J 



smt 



RA 



2(1 — cos 9) (sine + ismhq) 7 



JcJn 



(29) 
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Similar to Ref. |9|, it splits into factors depending either 
on massless or on massive coordinates. 

The ingredients of the action ((TJJ take a simple form: 

str(T 3 Q) 2 = 4(cos2c- cosh2g), 
str AQ = -4 [P + (1 - \)R^p] , 

where A = cos 8, P = cosh q — A cos c, R = cosh q — cos c, 
while the preexponent in (|13|) is the most involved: 



variables: 



str(fcAQ) = 4 [(cosh q + A cos c) + 2P 

+ (l + A)E£p + 2(l-A)E W £p + ---L 

where the omitted terms do not contribute to the average 
DOS. 

Let us expand the integrand in Eq. (I13[) in Grassmann 



F20 ip + F 02 pr\ + F 2 2 £ppv 



00 



The parametrization (|27f is singular at 9 = (Uc = 1) 
and at c = q = (f/ ro = 1). Therefore the integral (fT3|) 
will contain not only the regular contribution from the 
term F22, but also the contributions from the terms Fqq, 
F20 and F 2 which are finite due to the Parisi-Sourlas- 



Efetov-Wegner theorem [2l|, in complete analogy with 
the calculation of Ref. 0]. In particular, the terms with 
F o and F 2 o reproduce the C-class result (|T8|) . while the 
other two terms are responsible for the crossover. After 
some algebra, the general expression for the average DOS 
in the class CI - class C crossover takes the form 



(0(E)) = 1 



sin 2x 
2x 



"V 2 sinh 2 q — sin 2 c 

dc cos c 5 2 — sm ( x cose) sin(x coshg) exp[a(cos 2c — cosh2g)]. (30) 

-tt/2 sinh q + sin c 

I 



The DOS given by Eq. (f30|) is shown in Fig. 0] for sev- 
eral values of the parameter a. 

In the limit a = 0, Eq. (1301) reproduces the known 
result for class CI [6|: 

<f?ci(£)) = \ {x\3l(x) + Jl(x)] - J (x)Jx(x)) . (31) 

A small magnetic field (a < 1) changes the linear CI- 
class behavior (g(E)) oc E to the quadratic C-class be- 



havior (g(E)) oc E 2 at x ~ y/a: 



(g(E)) ^-xeri 



/8a 



(32) 



see inset in Fig. [4] A similar modification of the level 
repulsion exponent takes place in the orthogonal-unitary 
crossover, Eq. (p}. 

In the limit of large a, one finds 




FIG. 4: The average DOS, (g(E)}, in the class CI - class C 
crossover. The curves correspond to different values of the 
symmetry-breaking parameter: a — (class CI, tiny oscil- 
lations), 0.2, 1, and 00 (class C, large oscillations). Inset: 
(q(E)) for small deviations from class CI: a — 0, 0.01, 0.025, 
and 0.05. 



(0(E)) « 1 



sin 2x 



4a 



2x ^16a 2 + x 2 



(33) 



which reduces to the C-class result (1181) at a — > 00 



Conclusion. The purpose of this Letter was to study 
crossovers between superconducting symmetry classes. 
We have shown that the average DOS in the crossover re- 
gion between normal/superconducting symmetry classes 
with spin-rotation symmetry (Fig. [5]) can be calculated 
using Efetov's supersymmetric u-model of orthogonal 
symmetry. We have obtained exact expressions (I17[) and 
(|3H| for the DOS in the crossover regions between the 
classes A-C and CI-C. 
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